Astronomy & Astrophysics manuscript no. 6552 


© ESO 2008 


February 4, 2008 





The extra-mixing efficiency in very low metallicity RGB stars^ 



o 
o 

> 
O 



A. Recio-Blanco and P. de Lavemy 

Dpt. Cassiopee, Observatoire de la Cote d' Azur, CNRS/UMR 6202, BP4229, 06304 Nice cedex4, France 
e-mail: [arecio ; laverny] (Jobs-nice . fr 

Received 2006; accepted 

ABSTRACT 



> 



o 

6 



Aims. After the first dredge-up, low-mass Red Giant Branch (RGB) stars experience an extra-mixing episode that strongly affects the chemical 
abundances on their surface. This mixing occurs at the bump in the luminosity function. In this Letter we describe the efficiency of the 
extra-mixing in RGB stars found in very metal-poor globular clusters (GC). 

Methods. The VLT/ISAAC spectra of twenty stars located between the bump and the tip of the RGB in four GCs with metallicities between 
[Fe/H]=-1.2 and -2.5 dex were collected. The carbon isotopic ratios on their surface were derived from the second overtone (Av=2) bands of 
the CO molecule at 2.3/jm with the spectral synthesis method. 

Results. It is found that the carbon isotopic ratios of very metal-poor GC stars always reach the equilibrium value of the CNO cycle almost 
immediately above the bump in the luminosity function. No additional mixing episode at brighter luminosities and no variations with the 
clusters' metallicity were detected. The extra-mixing is therefore found to be very efficient in metal-poor low-mass RGB stars, in very good 
agreement with theoretical expectations. 
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1. Introduction 

During their ascent of the Red Giant Branch (RGB), low-mass 
stars undergo mixing processes that bring freshly synthetised 
nuclides to their surface. In standard stellar evolution theory 
(Iben. [T964l l. chemical changes in the surface are only expected 
to be caused by convective dilution during the first dredge-up 
(IDUP) occurring at the base of the RGB: the convective enve- 
lope deepens and enters regions akeady processed by the cen- 
tral H-burning phase, thereby altering the surface abundances. 
This process leads to a decrease in the carbon isotopic ra- 
tio ('^C/'^C) from the main sequence value (~90 in the solar 
case) to post IDUP values around 20-25 (Charbonnel 2003 or 
Denissenkov & Herwig 2004). In addition, the carbon abun- 
dance drops, while the nitrogen one increases. Oxygen and 
all heavier element abundances, however, would remain un- 
changed. According to the standard scenario, the surface abun- 
dances after the IDUP remain unaltered, as the convective 
envelope slowly withdraws outward in mass during the final 
stages of the RGB evolution. 

This canonical picture is, however, challenged by an in- 
creasing amount of observational data demonstrating its lim- 
ited validity; see, for instance. Brown & Wallerstein il989i and. 
more recently. Smith et al. 120001 and Gratton et al. I2000I I. In 



fact, in most of the field and globular cluster (GC) low-mass 
evolved stars, the observed conversion of '^C to '^C and '^N 
greatly exceeds the levels expected from standard stellar mod- 
els. More recently, Shetrone (120031) derived '^C/'^C ratios in 
low-metallicity GC RGB stars ([Fe/H] > -1.2) and found very 
low values, almost reaching the near-equilibrium value of the 
CNO cycle ('^C/'^C~ 3.5). This very clear evidence requires a 
non canonical mixing (named extra-mixing) between the shell 
and the bottom of the convective envelope (see for instance 
Charbonnel 120031 and Weiss l2006b . Most of the observations 
supporting this evolutionary scenario indicate that extra-mixing 
begins at the RGB luminosity-function bump. This bump ap- 
pears when the narrow burning H-shell reaches the sharp chem- 
ical discontinuity in the H-distribution profile caused by the 
deep penetration of the convective envelope. A drop in the stel- 
lar luminosity then occurs, revealed by a bump in the RGB lu- 
minosity distribution. This bump, predicted by Iben (I1968I I and 
detected by King et al. (1985) , occurs at brighter luminosities 
for more metal-poor GCs. 



* Based on observations collected with the VLT/ISAAC instrument 
at Paranal Observatory ESO (Chile) - Programme 75.D-0228A 



As already proposed by Charbonnel ( |1994|[T995] l, observa- 
tions strongly suggest that prior to the bump, the mean molec- 
ular weight gradient created during the IDUP acts as a bar- 
rier to any mixing below the convective envelope. After the 
bump, the gradient of molecular weight above the H-burning 
shell is much lower, allowing the extra-mixing to act. Now, the 
extra-mixing scenario is well admitted, but its exact physical 
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nature and efficiency are still objects of debate (see, for in- 
stance, Palacios et al. |20061 l. Furthermore, the extra-mixing ef- 
ficiency at very low metal content is still poorly known, as is its 
metallicity dependence. We therefore present here carbon iso- 
topic ratios of RGB stars in GC ten times metal-poorer (down 
to [Fe/H] ~ -2.5) than presented in Shetrone ( I2003I I. They were 
derived from low-resolution spectra around the second over- 
tone bands of the CO molecule at 2.3//m. The selection of the 
targets, the observations, and their reduction are presented in 
Sect. 2. Sect. 3 is devoted to the derivation of the '^C/'^'C ratios 
and, we discuss the extra-mixing efficiency in such metal-poor 
stars in Sect. 4. 

2. Target selection and observations 

We selected three GCs with [Fe/H] ratios about ten times lower 
than in previous works : NGC 6397, M30 (NGC 7099), and 
M15 (NGC 7078) with [Fe/H]=-2.1, -2.3, and -2.45 dex, re- 
spectively. For the purpose of a comparison with Shetrone 
(|2003] |. we also selected M4 (NGC 6121 with [Fe/H]=- 
1.2 dex). For all these clusters (see Table 1), we adopted the 
metallicity scale of Kraft & Ivans (2003). The observed K^- 
magnitudes of the Ml 5 and M4 bumps were taken from Cho 
& Lee 62002: FeiTaro et al '2000' and Valenti et al. "20041 re- 
port a consistent value for M15). For M30 and NGC 6397, 
we estimated their /Tj-bump magnitude from the relation de- 
rived by Valenti et al. (.2004) . from the photometric systems 
transformation of Carpenter (2001), and from distance moduli 
and extinctions in the Harris catalogue ( 119961 see next section 
for more details). This relation was checked for a very good 
agreement with the observed /T^-bump magnitudes of M15 and 
M4. The adopted /iT, -magnitudes of the GC bump are reported 
in Table 1. Typical uncertainties for M^-bump magnitudes are 
+0.1mag. (Valenti et al. I2004I I. Taking into account eiTors on 
distance moduli and extinctions, the total error for the calcu- 
lated /Tj-bump magnitudes is of the order of 0.15mag. It is 
around 0. Imag. for the observed ones (Cho & Lee l2002] l. 

Target selection was based on the photometry by Rosenberg 
et al. (I2000al ) for M4 and NGC 6397, Rosenberg et al.( l2000bt 
for M15, and Momany et al. ( I2004I I for M30. The colour- 
magnitude diagrams of each GC were used to select stars lo- 
cated between the bump (or slightly less luminous) up to the tip 
of the RGB. The observed stars are listed in Table 1 with the 
naming convention of their coiTesponding photometry. Some 
fainter targets of the sample were finally not included in this 
work due to an unclear identification of their CO bands. We 
also checked the radial velocity of each target, derived from the 
collected spectra, in order to confirm their GC membership. 

The observations were collected during the nights July 
25-27, 2005, with the ISAAC instrument on the VLT/ANTU 
(ESO, Chile). We used a slit width of 1", leading to a spec- 
tral resolution of about 3 000 over the domain 2.275-2.39/im. 
Objects were observed at two positions along the slit to remove 
the sky contribution (nodding technique). In most of the cases, 
the slit orientation was chosen in order to observe two objects 
simultaneously. Their position on the slit and the nodding/jitter 
parameters were defined to avoid spectra overlap. The number 
of nodding cycles varied from 3 to 5 and the detector integra- 



tion time from 3.5 to 300 sec, depending on the star magni- 
tude. The total exposure times (without overheads) ranged be- 
tween one minute for the brightest stars at the tip of M4 and 
NGC 6397, down to one hour for the faintest targets above the 
bumpofM30andM15. 

To remove the lines of the Earth's atmosphere, several tel- 
luric standards were observed at the same airmass and with the 
same instrumental configuration as the science targets, soon 
before or after them. The selected standards were featureless 
hot OB spectral type stars. The spectra were reduced with the 
ISAAC pipeline offered by ESO, and standard IRAF proce- 
dures were used for spectra extraction. Each target spectrum 
was then divided by the spectrum of its coiTesponding telluiic 
standard. In practice, several telluric spectra were tried in or- 
der to remove the telluric lines as fully as possible. Finally, 
the observed spectra were normalised as proposed by Shetrone 
(120031) . We first compute a synthetic spectrum with parameters 
and chemical abundances (see Sect. 3) as close as possible to 
the observed one. Then, the residual of the division of the ob- 
served spectrum by the synthetic one is fitted by a low-order 
spline, which is then used to normalise the observed spectrum. 
This procedure was performed iteratively until the best syn- 
thetic spectrum was found. 

3. Determination of the carbon isotopic ratios 

The carbon isotopic ratios were derived with the spectral syn- 
thesis method. Theoretical spectra were computed with the 
turbospectrum code (Alvarez & Plez 119981 and further im- 
provements by Plez). Model atmospheres were interpolated 
in a grid of new-generation MARCS models (Gustafsson et 
al. 120021 12006). These models are in LTE, spherical geom- 
etry and with an enhancement of the a-elements typical of 
such low-metallic environments ([ff/Fe] = -t-0.4). The CNO 
abundances of the Sun revised by Asplund et al. (I2005I I were 
adopted, and Solar abundances of other chemical species are 
from Grevesse & Sauval ( |1998l l. The adopted line list con- 
sists in atomic lines from the VALD database (Kupka et al., 
119991 ) and in CO lines, including its different isotopes, from 
(Goorvitch & Chackerian, 119941 ). Other molecular lines were 
found to be invisible in the studied domain at such low metallic- 
ities. Finally, synthetic spectra were broadened by convolution 
with a Gaussian to match the observed line widths. This line 
list was checked by fitting the high-resolution IR spectrum of 
Arcturus (Hinkle et al.. 119951) with stellar parameters and abun- 
dances from Peterson et al. ( fT993T l. We found '^c/'^C = 6+1 
for the Arcturus spectrum degraded to our adopted ISAAC res- 
olution, in very good agreement with previous determinations 
published since Griffin (119741 ). 

For each target star, the effective temperature was derived 
from the (V - K) and {J - K) colour indices and the Alonso 
scale (|1999ll200Tb . using the transformation equations between 
the different photometric systems from Carpenter (2001) and 
Bessel & Brett (1988). The IR magnitudes are from the 2MASS 
catalogue by Skrutskie et al. (20061). We estimated the extinc- 
tion in K from the GC foreground reddening reported by HaiTis 
d 1 9961 revised in 2003) and from the relations given in Cardelli 
et al. ( I1989I ). For M4, differential reddening was estimated from 
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Cudworth & Reed ( |1990| l but, in the effective temperature de- 
termination of its members, we favoured the estimate derived 
from the {J - K) colour, as it is less affected by the redden- 
ing. The surface gravity was calculated from the estimated ef- 
fective temperatures, distance moduli (Harris catalogue), and 
isochrones from Pietrinferni et al. ( I2004I I. Finally, we estimated 
the micro-turbulent velocity from the empirical relation for 
metal-poor giant stars derived by Pilachowski et al. ( 11996I) . 

Then, the carbon isotopic ratio of the target stars was de- 
rived following a procedure very close to the one adopted by 
Shetrone ( 120031 ) . The metal abundance of each star was as- 
sumed to be the same as the one of its parent cluster (in the 
Kraft & Ivans scale, [20031 see Table 1), and the oxygen abun- 
dance was adjusted from the a-enhancement of these GCs. We 
then varied the carbon abundance to match the ^^CO(4-2) band 
at 2.35/im. The '^C/'^C ratio was finally determined by fitting 
the '^CO(2-0) band at 2.345//m. We recall that the observed 
continuum was corrected when necessary with the new syn- 
thetic spectra (see Sect. 2) until a stable solution was achieved. 
We found that these two CO bands give the most reliable so- 
lution. It was indeed impossible to remove the telluric contam- 
ination accurately enough around the other '^CO band avail- 
able at 2.37/im. Finally, the adopted procedure was found to 
be insensitive to uncertainties in the stellar parameters, since 
both '^CO and '^^CO features vary with them in a similar way. 
Furthermore, because the '^CO band are not saturated in these 
very-low metallicity atmospheres, the uncertainty on the micro- 
turbulent velocity was found not to affect the results. Thus, 
once the '^CO band is fitted by assuming a given carbon abun- 
dance (and other stellar parameters and chemical abundances), 
only the '^C/'^C ratio can be varied to fit the '^CO feature. 

The derived carbon isotopic ratios are given in Table 1. 
The coiTesponding eiTors are completely dominated by (i) the 
signal-to-noise ratio of the observed stellar and telluric spectra, 
(ii) the removal of the telluric absorptions, and (iii) the nor- 
malisation of the observed spectrum (these last two points be- 
ing the dominant source of error). The error bars reported in 
Table 1 have been estimated by varying the '^C/'^C ratio, tak- 
ing all these uncertainties into account and, in particular, the 
continuum location around the '^CO and '^CO bands. For low 
'^C/'^C values, eiTors are found to be much smaller due to the 
greater strength of the '^^CO band. For higher '^C/'^C, this band 
can be very weak, leading to much larger error bars. 

4. Extra-mixing efficiency in very metal-poor stars 

Figure 1 presents our derived carbon isotopic ratios with re- 
spect to the evolutionary status of the targets above the lu- 
minosity function bump. Comparing the stellar /Tj-magnitude 
with the ATj-magnitude of the bump as in Fig. 1 has the ad- 
vantage of avoiding uncertainties in distance moduli and inter- 
stellar absorptions. Moreover, for the case of M4, this proce- 
dure diminishes any uncertainties due to diffential reddening. 
Errors in {Kg - ™'') are therefore dominated by the uncer- 
tainty on ^rf"™'' and are close to +0.15 mag in the worst cases 
(see Sect. 2). Of course, for a given cluster, the uncertainty on 
the relative positions of its members with respect to the bump 



Table 1. Target properties and derived carbon isotopic ratios. 
The adopted metallicities (in dex) of each GC are in the Kraft 
& Ivans scale (2003). See text for references on V, Ks, and 
A'f "'"''-magnitudes and discussion of their errors. 
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13.3 
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4+2 






665 


13.5 


2.3 


41 



luminosity is much fainter, since the 2MASS photometric error 
is about AKs - 0.04 mag. 

It can be seen that the carbon isotopic ratios are always 
found to be very low and close to or equal to the theoretical 
equilibrium value (within error bars), even immediately above 
the bump. In the case of NGC 6397 and M15, the '^c/i^c 
ratio is indeed found to always be smaller than 5 for a wide 
range of luminosities on the upper RGB, starting just above the 
bump. For M30, slightly higher values (around 6 in average) are 
found, but they are still compatible with the CNO equilibrium 
value and the extra-mixing scenario. We therefore do not detect 
any variation with metallicity: the extra-mixing appears to be 
very efficient whatever the value of [Fe/H] is, leading to a large 
decrease in the carbon isotopic ratio down to the equilibrium 
value. Unfortunately, we were unable to derive carbon isotopic 
ratios of stars fainter than the bump for these very metal-poor 
GCs in order to precisely define where the extra-mixing starts 
and what its exact duration is. This would require much larger 
SNR spectra with higher spectral resolution. Finally, there is a 
group of three M4 stars close to the bump but with rather low 
carbon isotopic ratios. These stars have therefore already expe- 
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Fig. 1. Evolution of the carbon isotopic ratio with respect to the stellar ^Tj-magnitude above the bump in the RGB luminosity 
function of each globular cluster. The adopted /Tj-bump magnitudes are reported in Table 1 . The maximum error on (Kg - /^f 
is shown in the upper right corner and, for a given cluster, the uncertainty on the relative stellar luminosities with respect to their 
bump is +0.04 mag. (see text for more details). The reported metallicity (in dex) of each GC is from Kraft & Ivans (I2003I I. For 
two stars (one in M4 and one in NGC 6397), we have only been able to derive upper limits for their carbon isotopic ratios. These 
limits are represented by a symbol located at the lower end of their error bar. 



rienced the extra-mixing episode, suggesting that the M4 bump 
could be slightly fainter (by about 0.1 mag.). 

On the other hand, our determinations for M4 stars can 
be compared to the work of Shetrone (12003I) . First of all, 
it can be seen that the carbon isotopic ratios in both stud- 
ies are very similar In addition, we observed two of the 
M4 RGB stars studied by Shetrone: the star M4#304 (named 
M4#4511 by Shetrone) and the star M4#456 (identified as 
M4#4507 by him). For the first one, we measure '^C/'^C=4+}, 
which is in very good agreement within the error bars 
with Shetrone's '^C/'^C=5;^2 determination. In contrary, we 
derived ^^C/^^C-l'^l for M4#456, whereas Shetrone found 
'^C/'-'C>15. We checked this discrepancy carefully and are 
fairly confident of our '^C/'-'C determination since (i) the SNR 
of our spectrum appears to be better than Shetrone's (see his 
Fig. 1) and (ii) M4#456 is slightly more luminous than the 
bump, so that it should have already experienced the extra- 
mixing episode as confirmed by our rather low estimate of its 



carbon isotopic ratio. Nevertheless, the same star is slightly 
fainter than Shetrone's M4 bump estimate. This disagreement 
could be explained by the uncertainties in the M4 distance 
modulus, which is needed to calculate the absolute magnitude 
values that he uses. Finally, Shetrone proposes that a small de- 
cline with increasing luminosity (from just above the bump up 
to the tip) was present in his data, suggesting a continued mix- 
ing. The rather small statistics per individual GC, together with 
the reported uncertainties on '^C/'^C, do not allow us to confi- 
dently explore the existence of this trend. 

Our conclusions can also be discussed with respect to the 
results of Gratton et al. (I2000I I who derived ^^C/'^C- 6^^ up to 
8^2 in nine field upper-RGB stars with metallicity lower than - 
2 dex. Although some of their estimates agree with ours, others 
seem to be higher than the CNO cycle equilibrium value. This 
discrepancy could result from the uncertainty on the precise lo- 
cation of the onset of the extra-mixing episode (i.e. the bump 
for GC stars) for field stars. Indeed, the maximum inward pene- 
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tration of the convective envelope occurs at higher luminosities 
for higher stellar masses or lower metallicity. In the Gratton 
et al. sample with stars with different masses and metallici- 
ties, this penetration therefore occurs at different luminosities. 
Thus, these RGB stars with slightly higher carbon isotopic ra- 
tios could still be experiencing their extra-mixing episode (in a 
GC, they would be found at or just above the bump). 

Finally, it has to be pointed out that Carretta et al. (I2005I I 
were the first to report determinations of the carbon isotopic 
ratio in GC subgiant stars. They surprisingly found rather low 
'^C/'^C ratios (mean value around 8) for these not very evolved 
stars that have still not experienced the IDUP. They convinc- 
ingly interpret their results by invoking the pollution of the 
analysed stars by previously evolved intermediate-mass AGB 
stars. If confirmed, this would reject the commonly admitted 
assumption that, before the bump, the carbon isotopic ratios in 
GC stars are still high. In consequence, the extra-mixing would 
not need to be as efficient as proposed by theoretical scenarios. 

On the other hand, our derived carbon isotopic ratios for 
low-mass very metal poor RGB stars are in very good agree- 
ment with extra-mixing estimates. For instance, Charbonnel 
(1995; see also Denissenkov & Vandenberg 12003 1) proposes 
that rotation induced extra-mixing leads to a decrease of 
i^C/i^C down to 3.5 for stars with 0.8 Mo and Z^IO"'* (i.e. 
close to the properties of our targets). In such low-metallicity 
stars, the mean molecular barrier, which inhibits the extra- 
mixing, would be much lower. That already leads to a very low 
'^C/'-'C ratio just above the bump luminosity as observed in our 
work. However, although the agreement appears very convinc- 
ing, these models are based on a parametric approach. More 
physically realistic simulations of rotating low-mass stars have 
been performed recently by Palacios et al. (2006) and show that 
self-consistent models, where the rotational transport is treated 
within the stellar evolution code, lead to chemical variations 
that are too small at the stellar surface. Open questions there- 
fore still remain in the modelling of rotating, evolved low-mass 
stars to reproduce the evolution of the observed chemical abun- 
dances on their surface. 

In summary, we therefore confirm that the extra-mixing 
scenario is efficient in any stellar system: field stars (Gratton 
et al., I2000I I. low (Shetrone, I2003I I. and very-low (this work) 
metallicity GC. Stars brighter than the RGB luminosity 
function bump already show the signatures of that mixing. 
Furthermore, within the uncertainties, the extra-mixing effi- 
ciency is found to be independent of the clusters metal content. 
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